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Supersonic turbulence in molecular clouds is a key agent in generating density enhancements that may subsequently go on to form 
stars. The stronger the turbulence - the higher the Mach number - the more extreme the density fluctuations are expected to be. 
Numerical models predict an increase in density variance, a 1 , , with rms Mach number, M of the form: a 1 , = b 2 M 2 , where b is 

r J PlPO PlPO 

a numerically-estimated parameter, and this prediction forms the basis of a large number of analytic models of star formation. We 
provide an estimate of the parameter b from 13 CO J= 1—0 spectral line imaging observations and extinction mapping of the Taurus 
molecular cloud, using a recently developed technique that needs information contained solely in the projected column density field 
to calculate a 2 , . When this is combined with a measurement of the rms Mach number, M, we are able to estimate b. We find 

PIPO ' ' 

b = 0.48+o [p which is consistent with typical numerical estimates, and is characteristic of turbulent driving that includes a mixture 
of solenoidal and compressive modes. More conservatively, we constrain b to lie in the range 0.3-0.8, depending on the influence of 
sub-resolution structure and the role of diffuse atomic material in the column density budget (accounting for sub-resolution variance 
results in higher values of b, while inclusion of more low column density material results in lower values of b\ the value b = 0.48 
applies to material which is predominantly molecular, with no correction for sub-resolution variance). We also report a break in the 
Taurus column density power spectrum at a scale of ~ 1 pc, and find that the break is associated with anisotropy in the power spectrum. 
The break is observed in both 13 CO and dust extinction power spectra, which, remarkably, are effectively identical despite detailed 
spatial differences between the 13 CO and dust extinction maps. 
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1. Introduction 

Q\ ' Recent years have seen a proliferation of analytical models of 
C*~) . star formation that provide prescriptions for star formation rates 
CN| ' and initial mass functions based on physical properties of molec- 
^ | ular clouds (Padoan & Nordlund 2002; Krumholz & McKee 
(Sj ■ 2005; Elmegreen 2008; Hennebelle & Chabrier 2008; Padoan 
O ' & Nordlund 2009; Hennebelle & Chabrier 2009). While these 
| models differ in their details, they are all fundamentally based 
t-H . on the same increasingly influential idea that has emerged from 
^ numerical models: that the density PDF is lognormal in form 
• i-h , for isothermal gas (Vazquez-Semadeni 1994), with the normal - 
■ ized density variance increasing with the rms Mach number 
H I ( cr p/ Po = b 2 M 2 where p is the density, po is the mean density, M 
- - ■ is the 3D rms Mach number, and b is a numerically determined 
parameter; Padoan, Nordlund, & Jones 1997). 

There is some uncertainty on the value of b. Padoan et 
al (1997) propose b - 0.5 in 3D, while Passot & Vazquez- 
Semadeni (1998) found b — 1 using ID simulations. Federrath, 
Klessen, & Schmidt (2008) suggested that b = 1 /3 for solenoidal 
(divergence-free) forcing and b = 1 for compressive (curl-free) 
forcing in 3D. A value ofb = 0.25 was recently found by Rritsuk 
et al (2007) in numerical simulations that employed a mixture of 
compressive and solenoidal forcing. Lemaster & Stone (2008) 
found an non-linear relation between <x 2 , and M 2 , but it is very 

PIPO J 

similar to a linear relation with b — 1/3 over the range of Mach 



numbers they analyzed. They also found that magnetic fields ap- 
pear to have only a weak effect on the o 3- . g - M 2 relation, and 
noted that the relation may be different in conditions of decaying 
turbulence. 
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Currently, observational information on the value of b, or 
the linearity of the - M 2 relation, is very sparse. A major 
obstacle is the inaccessibility of the 3D density field. Recently, 
Brunt, Federrath, & Price (2009; BFP) have developed and tested 
a method of calculating cr 2 p i po from information contained solely 
in the projected column density field, which we apply to the 
Taurus molecular cloud in this paper. Using a technique similar 
to that of BFP, Padoan, Jones, & Nordlund (1997) have previ- 
ously estimated a value of b — 0.5 (M =10, o- p / Po = 5) for the 
IC 5146 molecular cloud. In sub-regions of the Perseus molecu- 
lar cloud Goodman, Pineda, & Schnee (2009) found no obvious 
relation between Mach number and normalized column density 
variance, cr^m - It has been suggested by Federrath et al (2008) 
that this could be due to differing levels of compressive forcing, 
but it may also be due to differing proportions of o"p/ Po projected 
into cr 2 N i No (BFP). Alternatively, this may imply that there is no 
obvious relation between a 2 , and M 2 . 

P/Po 

The aim of this paper is to constrain the cr 2 / po - M 2 relation 
observationally. In Section 2, we briefly describe the BFP tech- 
nique, and in Section 3 we apply it to spectral line imaging ob- 
servations (Narayanan et al 2008; Goldsmith et al 2008) and dust 
extinction mapping (Froebrich et al 2007) of the Taurus molec- 
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Fig. 1. Integrated intensity map of the 13 CO J=l-0 line over the 
velocity range [0, 12] km s _1 in the Taurus molecular cloud. 



ular cloud to establish an observational estimate of b. Section 4 
provides a summary. 

2. Measuring the 3D Density Variance 

We do not have access to the 3D density field, p(x,y,z) in a 
molecular cloud to measure o 2 j g directly, but must instead rely 
on information contained the projected 2D column density field, 
N(x,y). One can use Parseval's Theorem to relate the observed 
variance in the normalized column density field, cr^, N to the 
sum of its power spectrum, PN/r? (k x ,k y ). It can be shown that 
PN/N (k x ,k y ) is proportional to the k, = cut through the 3D 
power spectrum, P p / Po (k x , k y , k z ). Since cr 2 ^ ^ can in turn be re- 
lated, again through Parseval's Theorem, to the sum of its power 
spectrum, P p / Po , it is possible, assuming an isotropic density 
field, to calculate the ratio, R = cr 2 , /cr 2 , from measurements 

' N/No' pi pa 

made on the column density field alone. It was shown by BFP 
that R is given by: 



R = 



A/2 A/2 

I Z ( p N/ No )(k)B 2 (k)]-P N/No (0)B 2 (0) 

k x =-A/2+l k y =-A/2+\ 



A/2 A/2 A/2 

[ Z Z Z ( p N/ No )(k)B 2 (k)]-P N/No (0)B 2 (0) 

k x =-A/2+l k y =-A/2+l k : =-A/2+l 



where A is the number of pixels along each axis, (PN/N )(k) is 
the azimuthally-averaged power spectrum of N/No, and B 2 (k) is 
the square of the Fourier space representation of the telescope 
beam pattern. In equationQ] it should be noted that {PN/N ){k) is 
the power spectrum of the column density field in the absence of 
beam-smoothing and instrumental noise. 

An observational measurement of cr 2 ^ and R can then yield 

°p/po' ^° ca l cu l ate tne power spectrum in equation Q] zero- 
padding of the field may be necessary, to reduce edge discon- 
tinuities and to make the field square. If a field of size A x x A y is 



zero-padded to produce a square field of size A px x A py then c 2 ^ 
should be calculated via: 



= l(l + ((,r: 



+ \)l 1 2 -\)IR p )-\, 



(2) 



where r\ = (A px A py /A x A y )i, and R p is the 2D-to-3D variance ratio 
calculated from the power spectrum of the padded field using 
equationQ] This assumes that the line-of-sight extent (in pixels) 
of the density field is A z = (A x A y ) 5 . 

The BFP method assumes isotropy in the 3D density field, so 
equation ([]]) must be applied with some caution. Using magne- 
tohydrodynamic turbulence simulations, BFP show that the as- 
sumption of isotropy is valid if the turbulence is super- Alfvenic 
(Ma > 1) or, failing this, is strongly supersonic (M > 10). The 
BFP method can be used to derive o 2 , to around 10% accuracy 
if these criteria are met, while up to a factor of 2 uncertainty 
may be expected in the sub-Alfvenic, low sonic Mach number 
regime. It should also be noted that the variance calculated at fi- 
nite resolution is necessarily a lower limit to the true variance. 
While this problem can only be addressed by higher resolution 
observations, estimates of the expected shortfall in the variance 
can be made from available information. 



3. Application to the Taurus Molecular Cloud 

3.1. Measurement of b using 13 CO data only 

We now apply the BFP method to 13 CO (J= 1—0) spectral line 
imaging observations of the Taurus molecular cloud. Figure Q] 
shows the 13 CO emission integrated over the velocity range 
[0, 12] km s _1 . In the construction of this map, we have re- 
moved the contribution of the error beam to the observed intensi- 
ties and expressed the resulting intensities on the corrected main 
beam scale, Tmb.c (Bensch, Stutzki, & Heithausen 2001; Brunt 
et al in prep.; Mottram & Brunt in prep.) The map is 2048 pix- 
els x 1529 pixels across, corresponding to 28 pc x 21 pc at a 
distance of 140 pc (Elias 1978). 

Following the procedure described in Section 2, we first es- 
timate the variance in the normalized projected field. For this 
initial analysis, we use the entire field as represented in FigureQ] 
with no thresholding of the intensities. We assume that the 13 CO 
integrated intensity, /o, is linearly proportional to the column 
density, N. The advantages and disadvantages of this assump- 
tion are discussed below in Section 3.2, along with an alter- 
native estimate of b using extinction data to calculate the nor- 
malized column density variance. Taking /13 oc N, we calculate 



N/N 



Io/h 



where 7o,i3 is the mean intensity. The observed 
variance in the field is the sum of the signal variance, cr 2 , and 
We measure cr 2 

= 2.39 (units areall (K km s" 1 ) 2 ). With 



the noise variance, <t 2 
°trf« = - 42 ' giving <r 



+ 2.81 and 

-1\2\ 



a measured 7o,i3 = 1-06 K km s we then find crj^y = 2.25. 

The power spectrum of the integrated intensity field is now 
calculated. We use a square field of size 2048 pixels x 2048 pix- 
els in which the map is embedded, and compute the power spec- 
trum using a Fast Fourier Transform. Application of tapers to 
smoothly roll-off the field edges had an insignificant effect on the 
result (see Brunt & Mac Low 2004). We calculate the noise floor 
using the power spectrum of an integrated map made over signal- 
free channels in the data, and this is subtracted from the power 
spectrum of the field. The central portion of the power spec- 
trum is shown in Figure|2]to demonstrate the level of anisotropy 
present. The black and white contours (as appropriate to the 
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Fig. 2. Central portion of the power spectrum of the Taurus 
13 CO J=l-0 integrated intensity map. The black/white contours 
(smoothed for clarity) show levels of equal power, while the red 
circles represent an isotropic power spectrum for reference. 




K 

Fig. 3. Zoom in on the central portion of the power spectrum of 
Figure |3] The white contours (smoothed for clarity) show levels 
of equal power, while the red circles represent an isotropic power 
spectrum for reference. 



greyscale level) show levels of equal power, obtained from a 
smoothed version of the power spectrum for clarity, and the 
red circles are overlayed to show what would be expected for a 
fully isotropic power spectrum. While Taurus is often considered 
"elongated" or "filamentary", the power spectrum of 13 CO emis- 
sion is in fact reasonably isotropic over most spatial frequencies. 
There is some evidence of anisotropy at the larger spatial scales 
(lower k), and we further demonstrate this in Figure [3] which 
shows a zoom in on the smaller k range of the power spectrum. 
Figure[3]shows that there is a preferred range in k for anisotropy 
to be present, occurring at spatial frequencies near k = 20. 

Figure |4] shows the angular average of the power spectrum, 
PN/N (k), with the noise floor subtracted and the beam pattern 
divided out. The heavy lines in Figure |4] display fitted power 
law forms (Pn/n<,Q0 x kr a ) to two distinct fc-ranges in the power 
spectrum. The dotted lines continue the fitted forms beyond their 
respective ranges for comparison. For k < 21, the best fitting 
spectral slope is a = 2.1, while for k > 21 we find a = 3.1. 
It is notable that the spectral break occurs at the preferred scale 
for anisotropy to be distinguishable in the power spectrum. The 
wavelength corresponding to a spatial frequency of A; = 21 is 
~ 0.5 degrees, or ~ 1.2 pc. Previously, Blitz & Williams (1997), 
using 13 CO data but not using the power spectrum, identified a 
characteristic scale of ~ 0.25-0.5 pc at which the structure of an- 
tenna temperature histograms changes notably. Hartmann (2002) 
also noted a characteristic separation scale of ~ 0.25 pc in the 
distribution of young low mass stars in Taurus. These scales are 
comparable to about 1 /4 the wavelength of the spectral break in 
the power spectrum. The orientation and scale of the anisotropy 
suggests that it arises from the repeated filament structure of 
Taurus, as suggested by Hartmann (2002). 

Figure [2] and Figure [3] show that the assumption of isotropy 
holds reasonably well in 2D, but this is not necessarily a guar- 
antee of isotropy in 3D, which, nevertheless, we must assume. 
The presence of a break in the power spectrum can be accom- 
modated by equation (Q]i which makes no assumptions, other 



than isotropy, on the form of the power spectrum. Evaluating 
equation (Q3, we find that R p = 0.029, where the subscript p 
reflects the fact that padding was employed in the power spec- 
trum calculation. Evaluating equation [2] with cr 2 ™ = 2.25 and 

77 = (Ap X ApylA x A y )T- = 1.157, we find that cr 2 / o = 74.23. The 
line-of-sight extent of the field assumed for this calculation is 
V28 x21 pc = 24.25 pc. 

Given the high estimated sonic Mach number in Taurus (see 
below) it is unlikely that large scale strong anisotropies like those 
seen in sub-Alfvenic conditions will be present due to mag- 
netic fields (BFP). Gravitational collapse along magnetic field 
lines could in principle generate anisotropy, and indeed the small 
anisotropy identified in the power spectrum is oriented as ex- 
pected in this case (Heyer et al 1987). 

To estimate the rms Mach number in Taurus, we now con- 
duct an analysis of the velocity field. Using the I3 CO spectral 
line averaged over the field (Figure |5}, we measure the line- 
of-sight velocity dispersion, cr 2 = (1.02 + 0.06) (km s~') 2 by 
fitting a Gaussian profile. The line-of-sight velocity dispersion 
includes contributions from all internal motions including those 
at the largest spatial scales arising from large scale turbulence 
(Ossenkopf & Mac Low 2002; Brunt 2003; Brunt, Heyer, & Mac 
Low 2009). Assuming isotropy, this implies a 3D velocity dis- 
persion of cr 2 = 3cr 2 = (3.06 ±0.18) (km s" 1 ) 2 . Taking the 
Taurus molecular gas to be at a temperature of 10 K (Goldsmith 
et al 2008), with a mean molecular weight of 2.72 (Hildebrand 
1983), the rms sonic Mach number is M — 17.6 ± 1.8, where 
the quoted error estimate is from spectral line fitting and uncer- 
tainties in the kinetic temperature (Goldsmith et al (2008) report 
kinetic temperatures for the majority of Taurus between 6 and 
12 K; we have taken an uncertainty of + 2 K). 

Combining these results, our observational estimate of b 
therefore is b — cr p / Po /M = 0.49 + 0.06, where we have applied 
an uncertainty of 10% to the estimated cr p/po (BFP 2009). The 
uncertainty arising from the assumption of isotropy is difficult 
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Fig. 4. Angular average of the 13 CO J=l-0 integrated inten- 
sity power spectrum (histogram) corrected for the noise floor 
and beam pattern. The solid lines show the fitted power laws 
P(k) oc k~ 2A and P(k) oc k~ 3 \ for the ranges k < 21 and k > 21 
respectively. The dashed lines show the continuation of the fitted 
power laws for reference. 




Fig. 5. Global mean spectral line profile for the Taurus 13 CO 
data (histogram). The heavy solid line is a Gaussian fit, with 
dispersion cr 2 = (1.02 + 0.06) (km s" 1 ) 2 . 



to quantify, and further progress on constraining b will require 
the analysis of many more molecular clouds. Extension of this 
analysis to a larger sample of clouds will allow a test of whether 
o 2 1 is indeed proportional to M 2 and may also allow investiga- 
tion of whether b changes due to varying degrees of compressive 
forcing of the turbulence (Federrath et al 2008; Federrath et al 
2009). The above estimate of the uncertainty in b accounts for 
measurement errors on M and known uncertainties in the calcu- 
lation of R from BFP. The true uncertainty is rather larger than 
this, as we discuss in the following sections. 




Fig. 6. Comparison of the CO integrated intensity map with 
the dust extinction map of Froebrich et al (2007). 



3.2. Measurement ofb using 13 CO and extinction data 

We now consider the consequences of the assumption of linear 
proportionality between the 13 CO integrated intensity and the 
column density. Goodman et al (2009; see also Pineda et al 2008) 
have recently evaluated different methods of measuring column 
densities in the Perseus molecular clouds. Their findings may be 
simply summarised as follows: 13 CO J=l-0 integrated intensity, 
7i3, is linearly proportional to Ay estimated by dust extinction 
over a limited column density range (over about a decade or so 
in /13), but is depressed by saturation and/or depletion in the high 
column density regime; in the low column density regime, /13 
is again depressed by lowered 13 CO abundance and/or subther- 
mal excitation, and is insensitive to column densities below some 
threshold, as evidenced by an offset term in the /13-Ay relation. 
While all these factors impact on the point-to-point reliability of 
using 7) 3 to derive N, their effects on o 2 , ,, and therefore on 

the estimated <r 2 , are of more relevance here. 

To empirically investigate the differences between using dust 
extinction and I3 CO integrated intensity as a measure of column 
density, we re-derived o" 2 ^ using the dust extinction map gen- 
erated from 2MASS data by Froebrich et al (2007). In Figure 
we compare, on the same coordinate grid, the the dust extinction 
map with the 13 CO integrated intensity map convolved to the 
resolution of the extinction map (4 arcminutes). The greyscale 
limits of each field were chosen to represent the fitted relation 
between the 13 CO integrated intensity, /13, and the dust extinc- 
tion, Ay, described below. While the overall spatial distributions 
of 7; 3 and Ay are broadly similar, there are detailed differences 
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Fig. 7. Plot of dust extinction versus 13 CO integrated intensity, 
lis, for the fields in Figure [6] The solid lines are the fitted rela- 
tions to A v - A v ,d + CI o listed in Table [1] 

in places, and column density traced by Ay is present on the pe- 
riphery of the cloud which is not traced by /13. 

The quantitative relation between /13 and Ay is shown in 
Figure|7] To these data, we fitted a linear relationship of the form 
Ay — Ay j) + CI b, where Ayp and C are constants. We use three 
different regression methods (Isobe et al 1990) and the fitted pa- 
rameters, Ayj) and C, are listed in Table[T] and the fitted lines are 
overlayed in Figure|7] While the scatter around these relations is 
quite large there is only a small number of positions in which 
significant saturation in /13 is obvious. 

Using the measured Ay field to calculate cr 2 ^^, under the 
assumption that Ay is proportional to column density, N, and 
accounting for a noise variance of 0.04 (Froebrich et al 2007), 
we find that cr 2 ^ = 0.70. This is significantly lower than 

°~n/n ~ 2-25 estimated through the 13 CO integrated intensities. 
There are two contributing factors to the lower value obtained 
using Ay. The first is simply that the fields shown in Figure [6] 
are of lower angular resolution than the field shown in Figure [T] 
A calculation of the normalized column density variance using 
the lowered resolution /13 map results in cr^/Wo = 1 -98, which 
is lower than that calculated from the higher resolution field, but 
still a factor of 2.83 greater than that calculated from the Ay field 
at the same angular resolution. (Further discussion of the resolu- 
tion dependence is given below and in Section 3.3). 

The second, and dominant, factor contributing to the lower 
°~n/n b m e asure d from the Ay field is the presence of the offset 
term, Ayo in the /13 - Ay relations given above. For reference, 
note that the addition of a constant positive offset to N will raise 
A^o but leave cr 2 N unchanged, thereby lowering o" 2 ^ . The true 
situation is more complicated that this, since the low column 
density material is also structured (i.e. contributes variance) - 
and it is possible that the column density variance scales with 
column density (Lada et al 1994). It is questionable whether this 
material should be included in the column density budget. The 
excess column density (roughly quantified by Ayo) is associated 
with atomic gas and possibly also diffuse molecular gas (H2) 
with little or no I3 CO. 



Regression method 


Ayo 


(mag) 


Ctmag/OCkms- 1 )) 


A v on / 13 


0.579 


± 0.002 


0.565 ± 0.003 


/13 on Ay 


0.331 


± 0.004 


0.799 ± 0.004 


Bisector 


0.462 


± 0.003 


0.676 ± 0.003 



Material in this regime does not contribute to the 13 CO emis- 
sion, and therefore also not to the mean spectral line profile, from 
which the Mach number is estimated. If the material is predom- 
inantly atomic with a possible contribution from diffuse H2, it 
is likely to be physically warmer than the interior regions of the 
cloud traced by 13 CO, thereby complicating the assumption of 
isothermality employed in most of the numerical models and 
in our calculation of the Mach number. Similarly, I3 CO emis- 
sion from the subthermally excited (but also likely physically 
warmer) envelope regions will contribute to both to /13 and to 
the mean line profile, but at a reduced level in relation to its true 
column density. In addition, while the "molecular cloud" part of 
the extinction is confined to a small region of around 24 pc extent 
at a distance of 140 pc, how the "diffuse" part of the extinction 
is generated is rather less well defined, and may include contri- 
butions from all regions along the line-of-sight to Taurus (and 
beyond). Consequently, the extinction zero level is quite uncer- 
tain. We argue therefore that the selection of material by I3 CO 
is in fact beneficial to our aims: it preferentially selects the inte- 
rior regions of the molecular cloud, to which the assumption of 
isothermality is likely to apply, and it limits the extinction bud- 
get to a spatially-confined region at the distance of the cloud. 
Note also that material contributing to 7)3 (and therefore N, by 
assumption) contributes in the same proportion to the mean line 
profile. 

Inclusion of opacity corrections, using I2 CO data, also po- 
tentially cause as many problems as they solve. A spatially uni- 
form opacity factor obviously has no impact on the calculated 
°%in ^ ue to me normalization. In principle, the compression 
of /13 in the high column density regime can be alleviated, but 
this effect is likely to be small - particularly since the strongly- 
saturated column densities account for a small fraction of the 
data. Goodman et al (2009) discuss numerous problems in the 
I2 CO-based correction for opacity and excitation in the low col- 
umn density regime; we also point out here that this is further 
complicated by the fact that 12 CO and 13 CO trace different ma- 
terial because of different levels of self-shielding and radiative 
trapping. If opacity corrections to /13 are made, then for consis- 
tency, detailed (largely impractical) opacity corrections should 
be made to the line profiles from which the Mach number is cal- 
culated, although for this the point-of-diminishing-returns has 
certainly long since been reached. 

To investigate the effect of the inclusion of the diffuse com- 
ponent, we calculate tff/ No from the field A v - Ay - Ayj>. For 
the first of the fitted relations in Table [1] (Ayo = 0.579) we 
find cr 2 N i Na = 2.73; from the second (Ayo = 0.331) we find 
^N/Na ~ 1-36, and from the third, bisector method (Ayr, = 
0.462) we find o^/ N = 1.91. These values may be summarized 
in a combined measurement, with associated uncertainties, of 
^N/Na ~ 1-91-0 55' which is more compatible with <t^, No = 1.98 
calculated from the I3 CO field at the same angular resolution. 
This is not surprising, as A v is, to a good approximation, lin- 
early correlated with /13, and the value of the coefficient C is 
irrelevant as we use normalized variances. 
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Fig. 8. Power spectra of the fields in Figure [6] The light his- 
togram is the power spectrum of 13 CO integrated intensity; the 
heavy histogram is the power spectrum of dust extinction. Both 
fields were scaled to the same global variance of unity before 
calculating the power spectra. The power spectra were corrected 
for the beam pattern but not the noise floor, resulting in an up- 
turn at high wavenumber. The straight and dashed lines are the 
fits from Figure [4] 



The dominant contribution to the uncertainty therefore 
comes from the choice made for the appropriate treatment of 
the diffuse component. We argue that the use of 13 CO, or A y , 
to calculate o 2 .,,,, is better-motivated, but we cannot resolve this 
issue further at present. 

We have noted above that the cr\, iM calculated from 13 CO 
is, naturally, lower in the reduced resolution field. Ideally, one 
should use the highest resolution data available from which to 
calculate <J^// N , and this raises the question of how the Ay- 



calculated value of o%j N may behave at increased resolution. To 

investigate this, we calculate the the power spectra of the 13 CO 
and Ay fields of Figure |6l after scaling each to the same global 
variance of unity (the mean of the field is irrelevant). The an- 
gular averages of the power spectra are shown in Figure [8] In 
these spectra, we have applied "beam" corrections, but have not 
attempted removal of the noise floor. 

It is remarkable that, despite detailed differences between the 
two maps, the power spectra are almost indistinguishable. The 
fitted power-laws from Figure|4]are overlayed for reference - de- 
tailed correspondance between the two power spectra is notable 
all the way down to the noise floor(s) near k ss 100. Comparison 
of the low resolution 13 CO power spectrum to the fitted slope at 
high wavenumber from Figure Q] shows that the noise signature 
turns on very abruptly, and the power spectrum is very reliable 
up to this point - this is a useful benchmark to gauge the equiv- 
alence of the two low resolution spectra. Previously, Padoan et 
al (2006) have reported differences between 13 CO and extinc- 
tion power spectra in Taurus, which is not seen in our analysis. 
The origin of this discrepancy may lie in the different extinction 
maps used by Padoan et al (those of Cambresy 2002), but may 
also arise from the influence of the effective "beam" size that is 



not corrected for in their power spectra as the resolution is var- 
ied. 

The effective equivalence of the 13 CO and Ay power spec- 
tra support our suggestion above that it is in the approximately 
uniform "diffuse" component that the principal statistical dif- 
ferences (at least those relevant to our analysis) between 13 CO 
and Ay are manifest. This, in turn, suggests that the Ay-derived 
value of <rjji N that would be obtained at the higher resolution 

afforded by the 13 CO data will increase in a similar manner 
to that seen in the 13 CO - i.e. an increase from 1.98 to 2.25, 
or a factor of ~ 1.14. Applying this to the A v -based measure- 
ments, our estimate for cr^^ at the highest resolution available 



is crt 



2.18+!?? 3 ,. Converting this to the 3D variance using 



N/N " -i "-0.63 

equation [2] (with r\ 



1.157 and R 



0.029) we find that o 2 ^ 



= 72.1^jg-^. Taking into account the measured Mach number, 
M = 17.6 ± 1.8, and applying an additional 10% uncertainty 
on the 3D variance arising from the calculation of R, we arrive 
at a measured b = 0.48^g jj. (For reference, if is calculated 
from Ay without subtraction of the diffuse component, we find 
b = 0.32.) 



3.3. The effects of unresolved variance 

As we have noted above, the measured variance of a field ob- 
served at finite resolution is necessarily a lower limit to the true 
variance of the field. An obvious question then arises: what about 
structure (variance) at scales below the best available resolution? 

If we assume that the observed slope of the power spectrum 
at high wavenumbers continues to a maximum wavenumber, 
k c , beyond which no further structure is present, we can esti- 
mate the amount of variance not included in our calculation (see 
BFP). It has been suggested that k c may be determined by the 
sonic scale: turbulent motions are supersonic above this scale 
and subsonic below it (Vazquez-Semadeni, Ballesteros-Paredes, 
& Klessen 2003; Federrath et al 2009). Taking cr v cc L 05 
(Heyer & Brunt 2004), we estimate the sonic scale in Taurus 
as L s = Lo/M 2 * 0.08 pc where Lq * 25 pc is the linear size 
of the cloud from Figure 1, and M = 17.6 is the Mach num- 
ber from above. This is comparable to the spatial resolution of 
the data (0.03 pc), which suggests that our measurement of <x^ po 
does not significantly underestimate the true value if structure is 
suppressed below the sonic scale. Taking k c — > 00 on the other 
hand would result in a factor of ~ 2 underestimation of <x 2 , and 

P/Po 

therefore a factor of ~ V2 underestimation of b. Our value of b 
must therefore be considered as a lower limit: the power spec- 
trum (Figure [4]i is not well measured near the estimated sonic 
scale (k s ~ 350), and further high resolution investigations are 
needed to resolve this issue. 



4. Summary 

We have provided an estimate of the parameter b in the proposed 
b 2 M 2 relation for supersonic turbulence in isothermal 

gas, using 13 CO observations of the Taurus molecular cloud. 
Using 13 CO and 2MASS extinction data, we find b = 0.48+° Jf, 
which is consistent with the value originally proposed by Padoan 
et al (1997) and characteristic of turbulent forcing which in- 
cludes a mixture of both solenoidal and compressive modes 
(Federrath et al 2009). Our value of b is a lower limit if sig- 
nificant structure exists below the resolution of our observations 
(effectively, below the sonic scale). If diffuse material is included 
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in the column density budget then we find a somewhat lower 
value of b — 0.32, although in this case there are some ques- 
tions over the assumption of isothermality in the Mach number 
calculation. Our most conservative statement is therefore that b 
is constrained to lie in the range 0.3 < b < 0.8, which is com- 
parable to the range of current numerical estimates. However, 
the o 2 J - M 2 relation has been tested to only relatively low 
Mach numbers (M < 7) in comparison to that in Taurus, so 
further numerical exploration of the high Mach number regime 
should be carried out. In principle, gravitational amplification 
of turbulently-generated density enhancements could increase 
a 2 " i at roughly constant Mach number, with gravity possibly 
acting analogously to compressive forcing. Numerical investiga- 
tion of the role of gravity in the cr 2 ^ po -M 2 would be worthwhile. 

Further analysis of a larger sample of molecular clouds is 
needed to investigate the linearity of the o 1 ^ - M 2 relation 
and possible variations in b. As our main sources of uncertainty 
arise because of questions over how to treat the diffuse (likely 
atomic) regime and how to account for unresolved density vari- 
ance, further progress on improving the accuracy of measur- 
ing b ultimately must involve answering what is meant by "a 
cloud" (Ballesteros-Paredes, Vazquez-Semadeni, & Scalo 1999) 
and down to what spatial scale is significant structure likely to 
be present? 
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Appendix A: Zero Padding 

Brunt, Federrath, & Price (2009) considered an observable 2D 
field, Fi, which is produced by averaging a 3D field, F 3, over the 
line-of-sight axis. The field F2 was considered to be distributed 
in a square region with a scale ratio (number of pixels along each 
axis) of A. Assuming isotropy, the 3D field then is distributed in 
a cubical region, again of scale ratio A. 

A column density field, N, is produced by integrating the 
density field, p, along the line-of-sight, rather than averaging. 
To conform to the conditions set out by BFP, it is necessary 
to work with the normalised column density field and density 
field, which are obtained by dividing each field by their respec- 
tive mean values, /Vo and po. The variances calculated from these 
fields are <t^ /JVo and cr 1 j respectively. The latter is the quantity 

required to test the theoretical prediction: o- 2 p ^ pa = b 2 M 2 . 

Consider a 2D field, F2, with mean value F%,o, which is the 
projection of a 3D field, F3, with mean F3J0. The normalised 
variance, cr 2 , = 0% / F > can be calculated for a square image, of 
scale ratio A, via: 



References 

Ballesteros-Paredes, J., Vazquez-Semadeni, E., & Scalo, J. 1999, ApJ, 515, 286 
Bensch, F, Stutzki, J.. & Heithausen. A., 2001, A&A, 365, 285 
Blitz, L., & Willams, J. P., 1997, ApJ, 488, L145 
Brunt, C. M., 2003, ApJ, 583, 280 

Brunt, C. M., Heyer, M. H., & Mac Low, 2009, A&A, 504, 883 
Brunt, C. M., & Mac Low, 2004, ApJ, 604, 196 

Brunt, C. M„ F ederrath, C, & Price, D. J., 2010, MNRAS in press, 

arXiv:1001.1046 (BFP) 
Elias, J. Ft., 1978, ApJ, 224, 857 
Elmegreen, B. G, 2008, ApJ, 672, 1006 

Federrath, C, Klessen, R. S., & Schmidt, W., 2008, ApJ, 688, 79 
Federrath, C, Duval, J., Klessen, R., Schmidt, W., & Mac Low, M.-M., 2009, 
ai'Xiv, 0905.1060 

Froebrich, D., Murphy, G. C, Smith, M. D.. Walsh, J., & del Burgo, C, 2007, 

MNRAS, 378, 1447 
Goldsmith, P. F, Heyer, M. H., Narayanan, G, Snell, R. L., Li, D., & Brunt, 

C. M., 2008, ApJ, 680, 428 
Goodman, A. A., Pineda, J. E., & Schnee, S. L., 2009, ApJ, 692, 91 
Hartmann, L., 2002, ApJ, 578, 914 
Hennebelle, P., & Chabrier, G, 2008, ApJ, 684, 395 
Hennebelle, P., & Chabrier, G, 2009, ApJ, 702, 1428 
Heyer, M. H. & Brunt, C. M., 2004, ApJ, 615, L45 
Heyer, M. H., Gong, H., Ostriker, E., & Brunt, C, 2008, ApJ, 680, 420 
Heyer, M. H., Vrba, F. J., Snell, R. L., Schloerb, F. P., Strom, S. E., Goldsmith, 

P. F, & Strom, K. M., 1987, ApJ, 321, 855 
Hildebrand, R. H., 1983, QJRAS, 24, 267 

Isobe. T., Feigelson, E. D., Akritas, M. G, & Babu, G. J., 1990, ApJ, 364, 104 

Kritsuk, A. G, Norman, M. L., Padoan, P., & Wagner, R., 2007, ApJ 665, 416 

Krumholz, M. R., & McKee, C. F, 2005, ApJ, 630, 250 

Lada, C. J., Lada, E. A., Clemens, D. P., & Bally, J., 1994, ApJ, 429, 694 

Larson, R. B., 1981, M.N.R.A.S, 194, 809 

Lemaster, M. N., & Stone, J. M., 2008, ApJ, 682, L97 



F 2 
1 2,0 



F 2 



(A.l) 



where: 
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and FziU j) is the value of F2 at the pixel (i, J). 

Similarly, the normalised variance of F3 is cr 2 = o~^? 3 / F30 , 
given by: 
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where: 
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F 3p ,o = (F 3p ) = " — ^ = (/^i = <F 3 >/77 3 .(A.14) 



Combining these results, we find: 
a% = {a\ + Drj 1 - 1 

o\ p = (p\ + \)t? -1. 



(A.15) 
(A. 16) 



Application of the BFP method to the zero-padded field 
yields the ratio of 2D-to-3D normalised variance: 



2 / 2 



(A. 17) 



Fig. A.l. Schematic illustration of zero-padding. 



F X o = (F 3 ) = 



ZZZ 

/=! 7=1 k=l 

A* 



where R p is calculated using the power spectrum of the zero- 
padded field, and we have identified this by the subscript p on 
R. The quantity of interest, however, is cr 2 which can be derived 
from the measured cr\ and R p via: 



(A.6) 



C 7| = -(l + ((cr? + l)^-l)//? p )-l, 



(A.18) 



and F 3 (i, j, k) is the value of F 3 at the pixel (z, j, k). 



obtained through combining equations l A. 1 51 |A 1 61 and lA.17l 
The BFP method assumes that the input image from which 



If F 2 is now zero-padded to a scale ratio of A p (shown ^ m pQwer spectmm is calculated is square. In the case that 
schematically in Figure [AD, the normalised variance, a\ p , of the obserV ed field is not square, but has dimensions A x x A y , zero- 



the resulting field, F 2p , is given by: 



CT 2p = 



F 2 

2p,0 



(Fl) 



padding to a square field of size A px x A py is required. If now r\ 
is defined as: 
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where now: 

An A n 



this allows cr 2 to be calculated via equation lA. 1 81 under the as- 
sumption that the line-of-sight extent of the field F 3 is: 
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A 2 , 
(F 2 )-5 = (F 2 )/rj 2 
A p 



(A.8) 
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A z — y/AxAy 

Obviously, it is desirable that A x 
isotropy to be best satisfied. 



(A.20) 

A v for the assumption of 



where we have defined r\ = A p /A, and noted that: 
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since F 2p is zero outside the region where F 2 is defined. 
Similarly: 
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